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ABSTRACT 

In this paper, we revisit the aquisition of angular momentum of galaxies by tidal shearing and 
compute the angular momentum variance cr^ as well as the angular momentum correlation 
function C^ir) from a peak-restricted Gaussian random process. This stochastic process de- 
scribing the initial conditions treats both the tidal shear as well as the inertia as dynamical 
fields and explicitly accounts for the discreteness of the inertia field. We describe the way in 
which the correlations in angular momentum result from an interplay of long-ranged corre- 
lations in the tidal shear, and short ranged correlations in the inertia field and which reflects 
the correlation between the eigensystems of these two symmetric tensors. We propose a new 
form of the angular momentum correlation function which is able to distinguish between par- 
allel and antiparallel alignment of angular momentum vectors, and comment on implications 
of intrinsic alignments for weak lensing measurements. We confirm the scaling L/M oc M^^^ 
and find the angular momentum distribution of Milky Way-sized haloes to be correlated on 
scales of ~ 1 Mpc//i. The correlation function can be well fitted by an empirical relation of 
the form C/.(r) oc exp(-[r/ro]''). 
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1 INTRODUCTION 



In the current paradigm, haloes acquire angula r momentum by tidal 
shear i ng from the ambien t matter di stribution jHeavens & Peacockl 
1988 1: 'Catelan & Theunj Il996al|bl : iLee & PeJBoOOl : iLed llOOft 



Lee & Park 2006; Schafer| | 2009l . for a rev i ew), w hich was first 
proposed by Hoylg a949h and ISciamj j 19551) . Tidal shear- 
well su pported by numerical simulations iWhit^ 198^ 



ing IS well supported by numencai simulations dWhit 
ISugerman et al.ll200d : ICatelan et al.ll200ll : iHahn et al.ll20g 



201C 



and leads to alignments of the angular momentum direction with 
the local tidal shear field. An important observational conse- 
quence of angular momentum alignments in the large-scale struc- 
ture, are induced intrinsic ellipticty alignm e nts betwee n neigh- 
bouring galaxies jCatelan & Porcianil 1200 ll : Ilin3 l2002l) . which 
can be expected to be a significant source of systematics 
in w e ak lensing surveys |Cro ft & Metzler 2000; Heavens et al. 



,2000; Crittenden et al." '200ll; Iffirata & SeliakI l2004i ; iKind \200i . 



Semboloni et al..,200&l and even galaxy surveys as they introduce 



selection efi^ects due to correlated angles of inclination of the galac- 
tic disks jRrause & Hiratalbol ih . By now, there is reliable obser- 
yational evidence o f tidal- s hearing induced ell i pticity correlations 
jMandelbaum et al.l l2006l : iHirata et alj |2007| : Ijones et al.l I2OI0I : 



e-mail: spirou@ita.uni-heidelberg.de 



iBlazek et alj201 ll) in particular with SPSS-da ta, and confirmations 
of these alignments in numerical simulations llAragon-Calvo et al.l 
I2OO7I ; iBetancort-Riio & Truiilldl2009l : ISchneider & Bridlell2010h . 

In this paper, we revisit the acquisition of angular momentum 
of cosmological objects in linear theory and recompute the corre- 
lation function of angular momenta. We restrict ourselves to linear 
structure formation, using the Zel'dovich mapping for the descrip- 
tion of the tidal shearing mechanism. In this paper, we hope to im- 
prove previous works on this topic in these aspects: 

(i) We employ an improved functional form for the correlation 
function which is able to distinguish between parallel and antipar- 
allel alignments of angular momenta, and which may assume neg- 
ative values for antiparallel orientation of the angular momentum 
vectors. This requires that the angular momentum correlation func- 
tion can not be a mere quadratic form in the tidal shear and inertia 
tensor fields, which are the relevant quantities for angular momen- 
tum build-up, but needs to be antisymmetric. 

(ii) Both the inertia and tidal shear fields will be consistently 
computed from a correlated Gaussian random process, such that 
the fields have consistent phase relations. The angular momentum 
correlation will reflect the different correlation lengths of the inertia 
and tidal shear fields. 

(iii) Treating both fields as dynami cal quantiti e s imp roves 
on the parameterisation introduced by iLee & Penl i2000l) and 
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ICrittenden et alj ( 1200 ih for the average misalignment between the 
eigensystems of both tensors and reflects changes in average mis- 
ahgnment with increasing distance. 

(iv) We explicitly take account of the discrete nature of the iner- 
tia field as the random process restricted to galaxy formation sites in 
the large-scale structure has a different weighting of certain inertia- 
shear combinations compared to that of continuous fields, i.e. the 
angular momentum distribution is biased. 

The theory is developed in Sect.|2] where we outline the Gaus- 
sian model used for determining the angular momentum correla- 
tions in the large-scale structure and where we propose an im- 
proved form of the angular momentum correlation function. In ad- 
dition, we comment on the influence of dark energy cosmologies 
on the angular momentum acquisition. The results are presented 
in Sect. |3] where we compute the angular momentum correlation 
function along with the angular momentum variance, and investi- 
gate their mass-dependence, followed by a discussion in Sect. |4] 
where we summarise our main results and comment on the con- 
sequences of the improved angular momentum model on intrinsic 
ellipticity correlations. 

Throughout, the cosmological model assumed is a spatially 
flat ACDM cosmology with Gaussian adiabatic initial perturba- 
tions in the cold dark matter distribution. Choices for the relevant 
parameter values are: Cl„, = 0.25 with Sli, = 0.04, Q.^ = 0.75, 
Hq = 100/ikm s-'Mpc"' with h = 0.72, n, = 1 and erg = 0.9. 



2 FORMALISM 

This section describes the Gaussian model used for deriving the an- 
gular momentum correlations in the large-scale structure. Sect. 12. 11 
explains how haloes acquire rotation by tidal shearing and relates 
the angular momentum L to the inertia I^p and gravitational shear 
'i'afi in the Zel'dovich-approximation. In Sect. 12.21 we outline a 
model for deriving the correlations of shear and inertia from the 
fluctuation statistics of the density field, based on a joint multi- 
variate Gaussian probability density. The covariances take a par- 
ticularly simple form if expressed in spherical coordinates, as ex- 
plained in Sect. 12.31 and we elaborate on the shape of the correla- 
tion matrices in Sect. 12.41 The correlation function {La(x)L„i(x')) 
of the angular momenta is determined in Sect. 12.51 by integrating 
out the Gaussian probability density restricted to peaks in the den- 
sity field. We discuss a technical issue, namely the misalignment 
in the shear and inertia eigensystems in Sects. IZ6l and l2.7l Due to 
the high dimensionality of the integration, we employ a numerical 
Monte-Carlo integration scheme, as explained in Sect. 12.81 



2.1 Acquisition of angular momentum by tidal shearing 

iDoroshkevici] imdl and'Whit^ h984h suggested that the angular 
momentum of galaxies originates from tidal torquing between the 
protogalactic region and the surrounding matter distribution prior 
to collapse. Assuming a non-spherical shape of the protogalactic 
region, the angular momentum grows at first order and linearly in 
time in Einstein-de Sitter universes, whereas in spherical regions, 
the acquisition of angular momentum is only a second order effect 
du e to convective matter streams on the boundary surface, as shown 
bv |Peeblesl ( IT96^ . 

Quite generally, the angular momentum L of a rotating mass 
distribution p(r, t) contained in the physical volume V is given by: 



L(t) = f dV (r-f)x v(r, t)p(r, t), (1) 
Jv 

where v(r, t) is the (rotational) velocity of the fluid element with 
density p(r,t) = {p)(l + 6(r,t)) at position r around the cen- 
tre of gravity f. In perturbation theory, (5 <c 1 and the den- 
sity field can be approximated by assuming a constant density 
( p) = Q. 

niPcrii inside the protoga lacti c region. Fol lowing White! 
jl984l) . ICateian & Theunj jl996ah and lCrittenden et alj ( l200lh . we 
describe the growth of perturbations on an expanding background 
in Lagrangian perturbation theory: The trajectory of dark matter 
particles in comoving coordinates is given by the Zel'dovich ap- 
proximation ( Zel'dovic h 1970): 



x(q, t) = q- D+{t)^^(q) ^ x = -O+VT, 



(2) 



which relates the initial particle positions q to the positions x at 
time t. The particle velocity x follows from the Zel'dovich-relation 
by differentiation by the time-variable. The growth function Z)+(0 
describes the homogeneous time evolution of the displacement field 
T and contains the influence of the particular dark energy model. 
In the Lagrangian frame, the expression for the angular momentum 
becomes 



L = poo" I A^q (x - x)xx =^ poo' I d'q{q-q)xx, (3) 

Jvl Jvl 

where the integration volume is defined in comoving coordinates 
as well. Assuming that the gradient VfC^) of the displacement 
field ^{q) does not vary much across the Lagrangian volume Vl, 
a second-order Taylor expansion in the vicinity of the centre of 
gravity q is applicable: 

d,^(q) - d„^>{q) + J]{q - q)p^„p, (4) 

The expansion coefficient is the tidal shear ^a-y at the point q: 

'i'.y(q) = d^dy^q), (5) 

because the Zel'dovich displacement field 4^ is related to gravita- 
tional potential O and can be computed as the solution to Poisson's 
equation AT = 6 from the cosmological density field S. The gradi- 
ent d„^(q) of the Zel'dovich potential displaces the protogalactic 
object, which is neglected in the further derivation, as we only trace 
differential advection velocities responsible for inducing rotation. 
Identifying the tensor of second moments of the mass distribution 
of the protogalactic object as the inertia /^^r, 



Jv, 



Ir. 



■poa^ I d^qiq- 

Jv, 



one obtains the final expression of the angular momentum L„ 



(6) 



(7) 



It is convenient to rewrite the time dependence of £)+ in terms of 
the scale factor a by dZ)+/d? = aH(a)dD+/da, yielding: 

da 



3 dD+ sr^ 
L„ - H(a) — — £„py > //jo-^ct 



(8) 



The theory of angular momentum acquisition by tidal shearing has 
been extended to nonlinear stages by using second order pertur- 
bation theory ( Catelan & Theuns 199 6b) and to include effec ts of 
non-Gaussian initial perturbations JCatelan & Theunslll997h . but 
for reasons of analytical computability, we restrict our model of 
angular momenta to the linear regime of structure formation of a 
Gaussian random field. 
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Figure 1. The influence of the dark energy model on the time evolution of 
angular momenta: Q(a) as a function of scale factor a for ACDM (wo = - 1 
and Wa = 0, solid line), for quintessence (wq = -2/3 and Wa = 0, dashed 
line) and a dark energy model with variable equation of state (wq = -2/3 
and w„ = 1/3, dash-dotted line). 



Fig.[T]compares the time evolution of the angular momentum 
in dark energy cosmologies with SCDM. We define the ratio 

gDE(«) 



(9) 



^fSCDM(fl) 

with q{a) = H(a)dD+/da (with D+(a) normalised to unity to- 
day and we pa rameterise th e dark e nergy equation of state with 
jChevalUer & Po larski 2 00ll ; lLinder & Jenkins.2003.) 



w(a) = Wo + (1 - a)Wa. 



(10) 



In SCDM these formulae simplify to H(a) = Hoa^^^^, D+(a) = a 
and consequently ^scdm = Hqw'^^. Fig.[T]suggests that the spin-up 
of haloes in dark energy models is significantly slower compared to 
SCDM, and the choice of the equation of state affects the time evo- 
lution significantly. The growth function D+(a) and its derivative 
dD+/da follows numerically as a solution to the growth equation, 



dlnH\ d 



din a 



—D^ = —n„,(a)D+(a), 
da Za-^ 



(11) 



in which the dark energy model affects the scaling of the Hubble 
function H(a) and of the matter density parameter n,„(a). In spa- 
tially flat dark energy cosmologies, the Hubble function H(a) = 
dlna/df is given by 



Hi 



% + (l- 



n,„) exp 3 



1 r dlna' 



(I + w(a')) 



(12) 



with the dark energy equation of state w(a). The value w = -I 
corresponds to the cosmological constant A. 

2.2 Gaussian model of the angular momentum correlations 

The goal of this section is to derive the 2-point correlation function 
of the angular momenta of objects that form at peaks in the cosmic 
density field. The quantities needed are the the tidal shear ^„f}{x) as 
well as the inertia Iap(x) of a peak region, which both can be related 
to the density field itself and its second derivatives. In cosmology, 
fluctuations in the distribution of matter are described by the over- 
density 6{x), which is defined as the fractional perturbation in the 



density field p{x), 6{x) = (p{x) - {p))l(p), with the average density 
(p) = n,„Pcrit. These perturbations are conveniently decomposed in 
Fourier modes S(k): 



S{k): 



S(x) exp(-ikx). 



(13) 



Specifying the power spectrum P{k) suffices to describe the statis- 
tical properties of a homogeneous and isotropic Gaussian random 
field: 



{d{k)S(kr} = (Infdiik - k')P{k). 



(14) 



The power spectrum, for which we choose the ansatz P(k) oc k"' ■ 
T^(k), is normalised to exhibit a variance of erg = 0.9 on scales of 
R = S Mpc/h by the relation: 



dkk'^W^(kR)P(k), 



(15) 



with a Fourier-transformed spherical top-hat for the filter function 
W(y), i.e. W(y) = 3 [sin(y) - y cos(y)] /y^. A common parameteri- 
sation for the sh ape of the transfer fun ction T(q) for CDM models 
was proposed bv lBardeenetal.Ul986l) : 



T{q) 



ln(l +2.34?) 



X 1 + 3 



2.34^ 

.S9q + (le.lqf + {5A6qf + (6.71#1 



(16) 



where the wave-ve ctor k is given in units o f the shape parameter F, 
first introduced bv lEfstathiouetal.1 jl992h . A convenient parame- 
terisation of the value of F as a functio n of the matter d ensity Q.„, 
and the baryonic density Qi, is given bv lSugivamal l ll995l) : 



: q(k) : 



k/Mpc-^h 



with F = n„,/i exp -fi, 



1 + 



(17) 



The mass scale of the objects of interest is set by imposing a 
smoothing on high spatial frequencies, where the diameter R of the 
isotropic filter function SR(k) corresponds to the size of the objects 
at the onset of collapse. For numerical reasons, we use a spherically 
symmetric Gaussian for S^ik): 



P(k) P(k)Sl(k) with Sdk) - 



exp(-k^R'/2). 



(18) 



In order to predict the correlation of the angular momenta L of 
objects which form at peaks in the Gaussian density field, one needs 
to relate the density gradient da(x), the second derivatives Sa/jix) 
and the tidal field ^^^(jc) to the density field Six): 



S(x) 



r d'k 
——r6(k)expiikx), 

J (2ny 



d6(x) 

dXa 

d^Six) 



r d^k 

' -zr-^kaS{k) sxp(ikx), 



- f^k 

J (2;r)3'^" 



kp6{k) exp(iAjc), 



(19) 



(20) 



(21) 



dxadxp 

The tidal shear follows from the solution of the Poisson equation 
Al'(x) = 6{x) linking the Zel'dovich potential ^{x) to the density 
field S(x): 

d^'Vix) f d^k Kkp ^ 



r d^k 
J 



-S(k)e,xp(ikx). 



(22) 



An important consequence of eqns. i2U and I l22t will be the fact 
that the angular momentum correlation is determined by two mech- 
anisms with differing correlation length: a short range correlation 
of the peak shapes and hence the inertia, and a long range correla- 
tion mediated by the tidal shear. 
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The joint distribution of the amplitudes of the density field, its 
derivati ves and the tidal shea r follows from a Gaussian probability 
density teardeen et al.lll986h : 

p(v)dv = p.xp I — v'\/'v\dv, (23) 

(2;r)"/2VdrtV \ 2 j 

where the quantities of interest at the point x have been arranged 
in a 15-dimensional vector v, i.e., 3 values for the density gradi- 
ent 6a(x), 6 values for the second derivatives Sapix) of the density 
field (due to the interchangability of the second derivatives) and 6 
values for the tidal shear ^^^(x), which is symmetric under index 
exchange as well. The covariance matrix V follows from the outer 
product Vjj = (vjv'j), = 1 ... 15. This probability density can 
be extended to include the field values v' at a second point x', 



p(w)dw 



(271)-^ VdeTW 



exp|--M''W 'wldw. 



(24) 



where the 30-dimensional vector w = {v, v') combines the vectors 
V and v' at the two points x and x' under consideration and the 
covariance matrix W, W;^ = (WjW*), = 1 ... 30, is defined in 
complete analogy. 

A peculiarity worthwhile mentioning is the fact that the den- 
sity field 6(x) is degenerate with the trace ti'i'„p(x) of the tidal shear 
because of Poisson's equation AT = trf^^ = 6. For that reason, the 
density field will appear in the above outlined random process as a 
derived quantity, whereas the entries of the tidal shear matrix will 
be drawn from the Gaussian distribution, with the peak restriction 
in place. 

The inertia tensor !„/} of an object forming at a peak in the 
density field at position Xp is related to t he second derivatives Sgg o f 
the density field at that particular point jCatelan & Theunslll996ah : 
In the eigenframe of the mass tensor -Sa/s = —dadpS at the peak, the 
density field can be approximated by a parabolic density profile, 



6(x) 



J 3 



(25) 



where or = 1, 2, 3 are the eigenvalues of the mass tensor. If the 
boundary dT of the peak region T is defined by the isodensity sur- 
face (5r = and if the peak height is expressed in units of the vari- 
ance cTo, Six) = vcTo, the boundary surface is given in the parabolic 
approximation by an ellipsoid equation: 



1, 



(26) 



where the semi-axes of the ellipsoid are related to the eigenval- 
ues Ag. by 



2vcro 



(27) 



The volume T of ellipsoidal peak region bounded by isodensity 
contour (5r = in the parabolic approximation is then given by: 



471 

r = y A,A,,A„ 



(28) 



which would immediately yield an estimate for the mass of the ob- 
ject: 



At: 

M = zyoF = T/oy A.,A,A-. 



(29) 



The inertia tensor follows from the second moments of the mass 
distribution, restricted to the volume T of the peak region, and is di- 



agonal in the mass tensor eigenframe. Carrying out the integration 
yields: 



hp = ^ T diag{Al + AlAl +Al, A] +A^). 



(30) 



The evolution of the density field is assumed to be homogeneous to 
first order, rjo = PoAq = <p)a\ with <p) = n,„p„i,. 



2.3 Describing the correlations in spherical coordinates 

Following the exami 3le of iRegos & Szaiavl j 19951) and 
iHeavens & ShethI ( Il999l) . we express the correlations between 
the density field, its derivatives and the tidal shear in spherical 
coordinates. The two peaks under consideration are assumed to be 
positioned on the z-axis, symmetric about the origin, and separated 
by a distance r, i.e. they have the coordinates x = (0, 0, +z/2) and 
x' = (0, 0, -z/2). The correlations take a particularly simple shape 
in the basis given by the set of dimensionless complex variables 

, iM'i r A^k 

yl,M=^^4^ —^k'^'-"6(k)Ye^(k)exp{ikxl (31) 

with k = k/k as the direction of the wave- vector k. cr^. are the 
weighted moments of the (smoothed) matter spectrum P{k): 



(32) 



The transformation between the physical frame and the y^^^-fmms 
for the scalar density field S{x) is given by 

o-oj!!o(^) = Six). (33) 
For the vectorial density gradient Saix, they read: 
(Tiy%ix) = ^6-ix), (34) 
aiyl^ix) = -^|JJ2{6Ax) + i6Ax)). (35) 
The tensor da/jix) can be determined from the v"^^^ -coefficients by: 

a-2y°oix) = -^/5/4{s,Ax) + 6y,ix)-26^ix)), (36) 
o-2j5,(x) = - VT572 (<5,,(x) + i(5,,(x)) , (37) 
o-iyUx) = +^|^5JS {6_Ax)-Syyix) + 2iS,„.ix)), (38) 



0"2Voo(-«) 



+ (<5,,(X) + Syyix) + S,-AX)) . 



(39) 



The relation linking the tidal shear 'i'a^ to the y"^^ -coefficients can 
be derived in complete analogy to eqns. ( 136b through i39\ . The y^^- 
coefficients of the tidal shear tensor field differ mainly by a factor 
of crilcTa from those of the mass tensor field, apart from the trace 
of the tidal shear: 

o-o>-2o'(a:) = -^{^Ax) + 'Vyyix)-2^-^ix)), (40) 
0-0^2,' (X) = - VT572 (l',,(x) + i'i'„(x)) , (41) 
arayl^ix) = + ^[151% 1^ Ax) - yyix) + 2W „ix)) , (42) 

(43) 



o-oyaoix) 



+ i;V,Ax)+'Vyyix) + ^:Ax)), 



emphasising the difference in correlation length between the den- 
sity field and the potential. 

The }'"jj^(j:)-basis inherits its symmetry under complex conju- 
gation from the spherical harmonics Y(,„'- 



yiix)' = i-iryijx\ 



(44) 
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which will become important at the stage of inverting the relations 
given above. Similarly to the vector v containing the physical vari- 
ables, the >'",jj(x)-coefiicients can be arranged in a vector y by map- 
ping the 3 indices n, £ and m to a new index i. The physical w-frame 
and the frame of the j-values are related by a linear unitary trans- 
formation. For clarity, we abbreviate y = y(x) and y' = y(x'). 



2.4 Shape of the correlation matrices 

As d emonstrated in the formalism proposed by ^ Regos & Szalavl 
( Il995,) , the correlation matrices needed in the Gaussian probabil- 
ity densities (eqns.[23l and 124b assume a particularly simple shape 
in the frame given by the y"„, -coefficients and can be expressed an- 
alytically in terms of moments of the dark matter power spectrum. 
The correlation matrix Y in this frame is defined as the expectation 
value 



(45) 



of the products of the elements in the vector {y, y') and can be split 
into two 15x15 submatrices: the auto-correlation matrix A, defined 
as Ajj = <>',.v*> and the cross-correlation matrix C = C(r), given by 
dj = iyiy'j), which depends on the distance r between the two 
points X and x' , 

Y-( ^ 
\ A 



(46) 



where is the Hermitean adjoint of C. The transformation be- 
tween the v"^^j -frame and the physical frame is given by the complex 
matrix R acting on the vector v and resulting in the vector y, and 
by the matrix S, computing {y,y') from w = (w, v'). 



R 
R 



(47) 



The matrices R and S can be constructed from the relations be- 
tween the v"^jj(x)-coefficients and the physical variables 5{x), 6a{x), 
6afi(x) and 'i'afi(x) compiled in Sect. 12. 31 



2.4.1 Auto-correlation matrix 

The auto-correlation matrix A in the v"^^-frame is defined by: 

A=A';L,,' = <)'L(^)y;',„'(^)*>- (48) 

Inserting the Fourier expansion of the variables, replacing the vari- 
ance of the density field {5{k)5{k')') with the matter power spec- 
trum P{K), and using the orthogonality relation of the spherical 
harmonics Y(,„{k), 



(49) 



(50) 



J"df2 Kf,,„,(fc)Kf,„,(fc)* = 5f,f,(5„,|,„2, 
yields for the auto-correlation matrix: 

Ann A nn c c i ^ \n-n I -i-n+n r. r. 
U'mm' - ^rm"H'"mm' - {— ^ ) Oce'Omm', 

O" f+2iiO" t+ln' 

where the definition of the cr^-coefiicients in eqn. l |32| l was used for 
substituting the (2 j + 2)"" moments of the power spectrum P(k). The 
structure of the matrix A is remarkably simple: It is diagonal in the 
indices £ and m and the sign of its entirely real entries is determined 
by whether n - n' is an even or odd number. 



2.4.2 Cross-correlation matrix 

The cross-correlation matrix C(r) is defined analogously, 

C('-) = C™_, (r) = {yl„{x)y1'„,' (x'T)- (51) 

The steps in simplifying this expression consist in inserting the 
definition of the >'"jj,(j:)-coefficients, in replacing the variance 
{6(k)d(k'y) with the matter power spectrum P(k) and in expand- 
ing the Fourier wave exp(ikr), r = x - x',by virtue of Rayleigh's 
formula. 



sxp{ikr) = 4nY,i''jdkr) J] Yi^„ifyYi,„{k). 



(52) 



The integration over the three spherical harmonics Y(„Xk) can be 
simplified by inserting the definition o f the Wigner-3 j symbols 
( lMessiahlll962l : lAbramowitz et al.lll988h . 



dnYf,^,(k)Y(„„,(krYt,„„(k). 



(-1)" 



nl,(2£. + 1)1 



An 







£i 
nil 



£i £i 

-m2 ni^ 



(53) 



Further reduction is reached by taking advantage of the fact that 
both peaks are assumed to lie on the z-axis. 



YtMif) = Sm 



2L+ 1 
4n 



(54) 



which yields the final form of the cross-correlation matrix C(r): 

c+c 



(-1)- 



0-f+2nO-r 



/ '■■^-^ l)i^^' ' KLj+('+2(n+n'-i 



X V(2^ + 1)(2^' + 1) I 





£' 




(o 





o)( 



£' L 
—m 



(55) 



where the A:"'jV(^'')-weighted spectral moments are abbreviated 
with iff,,, (r), 



KUr) 



27f- / 



dkrj((kr)P(k). 



(56) 



jjjkr) are the spheric al Bessel functions of the first kind 
jAbramowitz etal.lll988l) . The C™^^^, -coefficients are always real: 
The Wigner 3 7-symbols are unequal to zero if L + £ - £' is even, 
in which case i^+^'-'' is a real number. Furthermore, the summation 
over L can be restricted to the range \ f - f'\ < L < £ + £' due 
to the triangle condition applied to the Wigner-3 j symbols. The 
cross-correlation matrix C can be brought to block diagonal shape 
by suitable arrangement of the y"^^^ -coefficients in the vector y, more 
specifically, by grouping coefficients with constant value of m and 
increasing the modulus of m with increasing index 

In contrast to the constant values in the matrix A, the entries of 
the matrix C depend on the distance r = \x - x'\ of the two points x 
and x' . The symmetry of the entries of C under interchange of the 
points X and x' is given by the relation 



(y'L(x)y"'nA^r} = (-D'"' W,„(^')X„,(x)). 



(57) 



Typical correlation coefficients C"",^^^^^^, as functions of separation 
r are depicted in Fig. |2] The smoothing scale R (c.f. eqn. |18t has 
been set to R = I Mpc/h and the density threshold was chosen as 
V = 2, in order to represent galaxies. With the choice of CI,,,, the 
smoothing of the power spectrum at scale R corresponds to a mass 
scale of Mscaic = YPcrit^«^' ^ 3.1 X 10"Mo//7. For illustration 
purposes, the covariance matrix C"",^^, has been transformed to 
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2 3 

distance r = 



4 6 6 7 

- x'l in units of R = 



1 Mpc/h 



Figure 2. Correlation coefficients C"",^^,{r) = (y'(Jx)y'f,^, (x')*> as a func- 
tion of distance r = |x - x'| in units of the filter scale R, transformed into 



a frame, where ffie auto-correlation matrix A" 
trix. The smoothing scale R was set to R 
galaxies. 



is equal to the unit ma- 
1 Mpc//i in order to represent 




1 2 3 4 5 6. 

distance r = |x - x'lin units of R ■ 



1 Mpc//i 



Figure 3. Con-elation coefficients C'-'!!, ,(>') = (v'l W/i, , (x')*> as a func- 
tion of distance r = |x - x'| in units of the filter scale R, transformed into a 
frame, where the auto-correlation matrix A™,^^^^^^, is diagonal. The smooth- 
ing scale R was set to = 1 Mpc/Zi in order to represent galaxies. 



a frame, where the auto-correlations A!]' 



are diagonal and nor- 



malised to unity. 

Fig. [3] shows the distance dependence of the entries of 
CTp, t{y) in a frame, where the auto-correlation matrix Af,, , is 

it' mm' ^ ' ' ll'mm' 

diagonal instead of equal to the unit matrix. Notable differences 
to C™,^^, in this frame compared to the frame where to auto- 
correlations are unity include the fact that C™,^^, indicates the elon- 
gation of the isoprobability contours, which assume semi-axis ra- 
tios of - 5 and a less oscillatory behaviour. 



2.5 Ansatz for the angular momentum correlations 

The correlation function of the angular momenta in the large-scale 
structure follows from integrating out the 30d Gaussian probability 



density p(M')dM', constraint to peak regions in the fluctuating density 
field. The regions to which the integration is restricted are required 
to exceed a threshold v in density, 5{x) > vctq, to be of vanishing 
density gradient, Sa(x) = 0, and of negative curvature, dapix) < 0. 

The number den sity of maxima in the density field is modeled 
by a point-proce ss teardeen et all 1 19861 ; iRegos & Szalavl 1 1991 
iHeavens&Shetril 19991) : 



"peak = 2 Sli(x - Xi). 



(58) 



Close to the maximum / at Xj a Taylor expansion of the density 
gradient Sa(x) is applicable: 

Saix) = 2 S„ii{Xi)(x - Xi)p- (59) 
With this expansion one obtains for the peak density: 

"peak = ^Sd i^apixdSa) = |det 6„ll\ 5o(<5„). (60) 

The relation takes a simpler shape when considering the eigensys- 
tem of the mass tensor -Sa/j: Being a symmetric tensor, it has the 
three real eigenvalues /!,,;= 1 ... 3 which allows to replace the 
determinant by the invariant quantity |/li/l2/i3|. 

The constraints can be combined in a mask C(i>), which is de- 
fined as a function on the vector v containing the derivates of the 
Gaussian random field under consideration: Peaks in the density 
field are defined as points with amplitudes in units of the variance 
o"o = exceeding a certain threshold v and exhibiting a vanish- 
ing gradient 6„ as well as negative curvature Sap' 

C(v) = 4 [S„(x)] lAiA.M Y] ®(^,) lS(x) - o-ov] . (61) 

0(x) denotes the Heaviside step-function. The peak density np^aki 
i.e. the expecation value for the number density of peaks in the 
fluctuating density field 6 which exceed a threshold vctq can then 
be derived from the multivariate Gaussian random process p(v)dv. 



_ f 



dv p{v) C(v), 



(62) 



which corresponds t o the integral of the differential peak density 
«peak(^)dv defined bv lBardeen et alj ( Il986l) . 

In analogy, the expectation value of the angular momentum 
restricted to peak regions in the density field can be obtained with: 



1 



dv p{v) C(v)L(v), 



(63) 



npe3k(> v) . 

and the variance of the angular momentum field can be computed 
using: 



f 

'Ipeak(> y) J 



dv p(v) C(vW(v). 



In both formulae, the normalisation factor n ' , 

' peak 



(64) 

accounts for the 



discreteness of the measured quantity. 

Generalisation of the above relations to include a second 
peak results in the correlation function {La(x)La'{x')) of the an- 
gular momenta, with the Gaussian probability density p(M')dM' = 
p(v, v')dvdv': 

{L„{x)LAx')) = (65) 
-^4 T r dwC(w) r dv'CW) L,{v)LA'v')p(v,v')- 

"peakO ^) -J 

In general, the thresholds y, v' imposed on the peaks are equal. 
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As derived in Sect. 12. II the angular momentum L„ depends on the 
product of the inertia tensor Ip^ and the tidal shear ^a-y. 



with (X^)' = +X^. In the derivation of the angular momentum L, 



(66) 



if the acquisition of angular momentum of a protogalactic object 
is described in the Zel'dovich approximation. For convenience, we 
introduce the matrix X with the components: 



(67) 



Then, the correlation of the angular momentum components be- 
comes: 



{La{x)La,(x')) = a^b\eal3ye„>p.y{Xpy{x)Xp,y{x')). 



(68) 



In the next step we replace the Id variance of the components L„ in 
the correlation function by the 3d variance of the full vector L by 
taking the trace of eqn. J68l >. 

Cdr) = tr{L„(x)LAx')) = {L„(x)L„(x')), (69) 

which has the advantage of being a coordinate-frame indepen- 
dent quantity and allows the usage of the relation E„pyE„p>y> = 
'i{Sfip>6yy> - SpySpiy"^ foT Teduciug the product of the two e„/jy- 
symbols: 

<L„(x)L„(x')) = 3a'Dl [(Xpyix')Xpy(x)} - {Xpy(x)Xyp(x')}] , (70) 

where the order of the indices in the last term is interchanged. In 
matrix notation, the correlation function C^ir) reads: 

Cdr) = tr{L(x)L'(x')) = Sa-^/J^tr [<X(x)X'(a:')) - <X(jc')X(x))] , (71) 

which is non- vanishing for general asymmetric matrices X^y due to 
the matrix transposition in the last term. Finally, one can express 
the angular momentum correlation arising from the Gaussian prob- 
ability density in the natural variable X, X^y = Ipa-'i'a-y'- 

tr<L(x)L'(x')) = 3^ — X 



(72) 



\^ AvC{v) J dv'C{v')tr[X(v)X'(v')-X{v')X(v)]p{v,v'). 

In the calculation outlined above we aim to avoid a decomposition 
of the tidal shear and inertia correlations according to 



(Kx)'¥ix) Kx'mx')) = P(\\'¥) P(\\V) CVix) '¥(x')). 



(73) 



which uses the difference in correlation lengths of the 1" and l-fields 
and is valid on scales on which the l-tensors are uncorrelated. 



2.6 Misalignment of the shear and inertia eigensystems 

The tensor X = l*P can be decomposed according to X = X^ + X" 
into an antisymmetric contribution X", defined via the commutator 
[1,1'], 

X-^i[l,«P], Xpy=^j]{lp^%,y-^>p^I^y), (74) 

(T 

with the symmetry (X")' = - <¥])' = i(TI - IT) = -X" under 
matrix transposition (I and *P are symmetric matrices) and into the 
corresponding symmetric matrix X^ by using the anticommutator 
{I, *?) between inertia I and tidal shear *P: 



X^ ^ i {I, , = 5 X (VT^, + ^,J.y) , 



(75) 



La — a~D+£„Ry y Ipcr^a-y — <^^L)-f-£„n 



(76) 



the permutation symbol e„^^ picks out the antisymmetric contribu- 
tion X", by virtue of £ai3y(Xpy + X^^) = X^^, because the contrac- 
tion of the symmetric tensor X^ with the antisymmetric permuta- 
tion symbol vanishes, e^pyX^^ = 0. 

Hence, the protogalactic objects will only acquire angular mo- 
mentum if the commutator X" between the inertia and the tidal 
shear is non-zero, which means that the inertia and shear tensors 
are not supposed to be simultaneously diagonisable, i.e. they are 
not allowed to have a common eigensystem. In order to capture this 
mechanism. iLee & PenI ( I2OO0I) and Crittenden et al. (2001) have 
used an effective, parameterised description of the average mis- 
alignment of the shear and inertia eigensystems, gauged with nu- 
merical w-body data. 

In the correlation function of the angular momenta Cdr) = 
tT{L'{x)L(x')) (c.f. eqn.l7U. the dependence of L on the commuta- 
tor X" translates into the asymmetric quadratic form (X(x)X(x') - 
X(jt:)X'(jc')) with the matrix transpose in the second term carrying 
the signal: A common eigensystem of the inertia and shear tensors, 
being both symmetric matrices, would have the consequence that 
X would be a symmetric matrix, X = X* = X', and the correlation 
function Cdx) would vanish. 

In contrast, the signal is maximised, if the shear and inertia 
eigensystems are unaligned, i.e. if X is purely antisymmetric, X = 
X" and X^ = 0. In that case X' = (X")' = -X" = -X and the angular 

(r) is truly quadratic: 



momentum auto-correlation function C""" 



(r) 



^4 ^2 n n 

■ 6—, — dvC(v) dv'C(v')ti\X(v)X'{v' 

nL.A> y)J J ^ 



d> v) , 
peak ^ ^ 

The corresponding 1 -point variance would be given by: 



(-^max) 



1ps^k{> v) J 



dv p(v)C(v)Ll,,(v), 



) .(77) 



(78) 



which illustrates the effect of partial alignment of the shear and 
inertia eigensystems, reducing the variance compared to the case 
where the shear and inertia eigensystems are maximally mis- 
aligned. 

The symmetric contribution X^, which measures the degree 
of alignment of the inertia and shear eigensystems, causes an 
anisotropic deformation of the protogalactic region during the 
course of linear structure formation prior to gravitational collapse. 
Consequently, the determination of ellipticity distributions is likely 
to be affected even in the stage of linear structure formation, and 
predic tions of triaxiality ba sed on peak shapes in Gaussian random 
fields teardeen et al.lll98^ could be refined using an adaptation of 
the formalism outlined above. 



2.7 Symmetry of the cancellation mechanism 

A possible objection concerning the symmetry of the correlation 
function Cdr), 



Cdr) oc (tr(X(jc)X'(x') - X(x')X(x))) , 



(79) 



might be that the mechanism comparing the alignments of the shear 
and inertia eigensystems outlined above is only present in the ma- 
trix X(x') at position x', which appears transposed in the second 
term, but not in the matrix X(x) at position x. The expression, 
however, can be reformulated using tr(X(jt:')X(jt:) - X(x)X(x')') = 
U(X(x')X(x) - (X(jt:)X(A:')')') = tr(X(x)X(x') - X{x')X{x)'), using 
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the properties of the trace tr(A + B) = tr(A) + tr(B), tr(AB) = tr(BA) 
and tr(A') = tr(A), as well as (A')' = A and (AB)' = B'A' such that 
the mechanism is present in the first tensor as well. 

Another objection might be that the correlation function Ctir) 
vanishes if one of the matrices X is symmetric. A symmetric shape 
of the matrix X, however, never occurs in the angular momentum 
build-up, because of the fact that symmetric matrices form a group 
under matrix multiplication: A symmetric matrix X could only have 
emerged from symmetric matrices I and *P with a common eigen- 
system due to a vanishing commutator [l,*P], making it impossible 
for the halo to acquire angular momentum as discussed earlier, and 
consequently, the correlation function has to be zero. 

2.8 Numerics of the constraint Gaussian integration 

The covariance matrix Y in the frame given by the y"^ -coefficients 
is transformed to the physical frame yielding the correlation ma- 
trix W. Then, we determine numerically the correlation functions 
Cl('') = ii{L'(x)L(x')) of the angular momentum L as a func- 
tion of distance r = \x - x'\, by carrying out the integration over 
the multivariate Gaussian probability density p(w)dw, subjected 
to being constraint to the peaks in the density field by the mask 
C(w) = C(w)C(w')- The inertia !„/} of an object forming at a peak 
in the density field is consistently derived from the local curvature 
dadp6{x) of the density field at the peak. 

Due to the high dimensionality the numerical constraint in- 
tegration is a difficult and time consuming task. This task is most 
efficiently addressed by exploiting our prior knowledge of the un- 
derlying probability density function: Since the distribution is just 
- tough very high dimensional ~ Gaussian, it is advantageous to 
sample the integral of eqn. M2\ directly instead of relying blindly 
on a common Monte-Carlo scheme. Generating samples which fol- 
low the distribution given in eqn. (|23} and eqn. ( I24t respectively is 
straightforward since this can be mapped onto the generation of 
unit Gaussian variates, for example via the Cholesky decomposi- 
tion of the corresponding covariance matrix. Gaussian (unit) vari- 
ates, however, can be obtained from a variety of very fast and effi- 
cient random number generators, e.g. the ziggurat method. 

Following the strategy of reducing the numerical integration 
of eqn. M2\ to a direct sampling process which only requires the 
generation of unit Gaussian variates is the most important sim- 
plification and acceleration in comparison to the use of standard 
Monte-Carlo techniques. In addition, our method allows to carry 
out the sampling in the physical frame where the constraints C can 
be evaluated most easily. Further acceleration can be achieved by 
using simple linear algebra in order to minimize the number of di- 
agonalisations of the mass tensor Map = dadp6. Instead of checking 
whether all eigenvalues of the mass tensor are negative we check its 
negative definiteness by combining the invariants tr(M"), « = 1,2 
and the determinant det(M). Finally, it is important to note that the 
actual dimensionality of the integration is reduced by the require- 
ment of being at a peak in the density field, i.e. 5o(V(5(x)). Thus, 
in our sampling process the derivatives of the density field are no 
degree of freedom but rather fixed (namely to be zero). This re- 
quires to adjust the overall normalisation of our constrained Gaus- 
sian sampling process by an appropriate factor, i.e. we have to use 
the normalisation of the unconstrained Gaussian distribution. 

In order to investigate the performance and reliability of our 
sampling strategy used to evaluate integrals like those of eqn. l l64t 
and eqn. Ml\ we computed in a first step the mean peak density 
given in eqn. ( 162) and compared our re sult with the analytical so- 
lution derived by ISardeen et alj ( Il986h . Increasing the number of 



-values with errors x 1 00 

> ,0.6469±0.001 




halo mass [Mq/Ii] 

Figure 4. Standard deviation a-^/u of the specific angular momentum as 
a function of halo mass, as obtained from numerical sampling (circles) 
including error bars (artificially enlarged by a factor of 100), along with 
a power-law fit to the data (solid line) and the theoretical expectation 
o"L/M M^^^ (dashed line). 



samples we could reproduce their analytical result to arbitrary accu- 
racy. We could also recover the limiting case v — > -cm. To provide 
an error estimate for our results of the angular momentum correla- 
tion functions given in the next section we carried out every sam- 
pling process several times seeding the random number generator 
differently and computed the resulting Itr error. For all numerical 
calculations we assured that the number of accepted samples is of 
the order of 10^ for each step in distance. 



3 RESULTS 

For visualisation, we define the specific angular momentum 



(80) 



normalised by the mass scale Mscaic = ^PaH^mR^ and divide out 
the Hubble-constant Hq defining the cosmological time scale such 
that the resulting quantity has units of a squared length scale and 
depends only on a single fundamental unit, for which we choose 
Mpc/A. 

As shown in lPeeblesI ( Il969h and lHeavens & Peacockl ( Il988l) . 
the standard deviation (Ti/m of this quantity scales oc M^'^, which is 
a valuable check for our numerical sampling code. In fact, as shown 
in Fig. |4l we recover this behaviour to a high degree of accuracy 
over two orders of magnitude in halo mass. A fit to the data points 
including their sampling error yields 



O" L/M 



O-L/M.O \ Mo I 



(81) 



with atiM = 11.3620-0 (Mpc/hy and Mo = 10^^ Mq/Ii, within 3% 
of the theoretical value. For convenience, we have divided L/M 
in Fig. |4]by ctq. Physically, the lo" value of L/M corresponds to 
a homogeneous sphere with the radius 1 Mpc/h and a mass of 
lO'^Mo//; revolving once every ~ 10' years. 

The correlation function Ci(r) of the specific angular momen- 
tum is depicted in Fig.|5]for a mass scale of IO'^Mq/A correspond- 
ing to a Milky Way-sized halo. It is decreasing rapidly after a dis- 
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Q correlation C^(r) 
fitting modei 



distance r [Mpc/h] 

Figure 5. Angular momentum correlation function Ci(r) for a Milky Way- 
sized halo of IO'^Mq//! (circles) including sampling errors, together with 
an empirical fitting formula of the form Ci^{r) oc exp(-[r / ro'f) (solid line). 



tance of ~ 1 Mpc/h corresponding to tiie correlation lengtii of 
the angular momentum field and assumes the asymptotic value of 
Ctir) — > cr^/ipcak for r — » in fulfilment of the Cauchy-Schwarz 
inequality with our values for the peak density. An empirical fit to 
the spectrum is given by 



Ci(r) = flexp(- [r/rof) 



(82) 



with a = (25010 ± 165)cr5 {Mpc/h)\ ro = (0.8628 ± 0.008) Mpc//i 
(both error bounds correspond to Icr) and p = 3/2. In addition, 
we find the angular momenta of haloes to be positively correlated, 
i.e. angular momentum vectors tend to be aligned in a parallel way, 
because the form of our correlation function would be able to dis- 
tinguish between parallel and antiparallel alignment. 



4 SUMMARY 

In this paper, we recompute the angular momentum correlation 
function arising from tidal torquing in an improved statistical 
model, which is based on a peak-restricted Gaussian random pro- 
cess. 

(i) Dark energy influences angular momentum aquisition by the 
time derivative H(a)dD+/da of the growth function D+(a). Angular 
momentum build-up in dark energy models is significantly slower 
compared to the SCDM-cosmology: At early times the angular mo- 
menta grow 30% slower in ACDM and 20% slower in classical 
quintessence with wq = -2/3, with qualitative differences in dark 
energy models with epoch-varying equation of state w(a). 

(ii) Due to the fact that the angular momentum arising in tidal 
shearing is proportional to the commutator X" = between 
the inertia I and the tidal shear *P the angular momentum field is 
sensitive to the relative misalignments of the principal axis systems 
of the tidal shear and inertia tensors. For that reason both fields are 
included in the Gaussian random process. 

(iii) The angular momentum correlation function and variance is 
computed from a peak-restricted Gaussian random process, due to 
the fact that galaxy formation is associated with local peaks in the 
density field. Technically, the peak restriction introduces a different 



weighting of the tidal shear-configurations compared to a continu- 
ous field. The inertia of a peak region is determined from a local 
parabolic density profile, and integrated within the boundary given 
by (5 = 0, which might be a too coarse approximation in the com- 
putation of the inertia, because of the weighting of mass elements 
with the square of the distance from the centre of gravity. 

(iv) The covariance matrices used for describing the relation be- 
tween the individual Gaussian derivates, are expressed in the y"^- 
basis (i.e. in spherical coordinates), and transformed to the phys- 
ical variables with a linear transformation. An important detail is 
the degeneracy of the density field 6 with the trace of the tidal 
shear 2t, d^da^ = A*? due to the Poisson equation. For that rea- 
son the density field itself, although it is used for deriving the rela- 
tions between all variables in the Gaussian random process, is not 
part of the random process itself, but derived from the tidal shear. 
Consequently the sampling can a priori not be restricted to values 
S > vctq, but 6 needs to be recomputed for each sampling point, 
and eventually be rejected, which results in a sampling efficiency 
amounting to ^ 2.5% for the 2-point function. 

(v) Computing the standard deviation from of the distribu- 
tion of the specific angular momentum L/M confirms the ex- 
pected scaling oc M^^^ prop osed by Peebles ( 1969) and derived by 
iHeavens & Peacock ( Il988f) in an analytical approach, and provides 
a valuable check for our code. 

(vi) The resulting correlation function Ci(r) of the specific an- 
gular momentum L = LKHoM^.^i^) (where the rescaling makes it 
possible to state the angular momentum in units of a single funda- 
mental unit) is determined by the interplay of two fields with dif- 
ferent correlation lengths: The inertia la/s exhibits short ranged cor- 
relations, and the tidal shear has the same correlation length as 
the density field. The resulting correlations in angular momentum 
have a range of about 1 Mpc/h for Milky Way-sized haloes, and 
can be fitted well with an empirical formula C/.(r) oc exp(-[r/ro/). 

Future investigations will include the application of the angu- 
lar momentum correlation function for computing ellipticity corre- 
lation functions, which play an important role in gravitational lens- 
ing: There, a common assumption are uncorrelated galaxy shapes, 
but coupled angular momenta give rise to the effect that the galac- 
tic disks of neighbouring galaxies are viewed under similar angles 
of inclination, such that their ellipticities are naturally correlated. 
The formalism outlined above can serve to compute the ellipticity 
correlation function, and furthermore, the intrinsic shape E- and B- 
mode correlation functions, Cf (f) and Cf{£) of the (tensor-valued) 
ellipticity field e - a paper about this topic is in preparation. 
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